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ABSTRACT 

A massive young star cluster, initially embedded in its parent molecular cloud, will spiral into the 

1/2 

Galactic Center from < 30m 6 pc during the life-time of its most massive stars, if the combined total 
mass is ~ 10 6 mg M Q . On its way inwards the system loses most of its mass to the strong tidal field, until 
the dense cluster core of high-mass stars is finally disrupted near the central black hole. A simple model 
is presented to argue that this scenario may under plausible conditions explain the observed location 
and rotation of the Galactic Center Hel stars. Accretion of star clusters into the Galactic Center could 
be recurrent, and play an important role in regulating the activity of Sgr A*. 

Subject headings: Galaxy: center - Galaxy: evolution - galaxies: nuclei - galaxies: star clusters - ISM: 
clouds - black hole physics 



1. INTRODUCTION 

The central parsec of the Galaxy contains a cluster of 
young stars, including some 15 very luminous Hel emis- 
sion line stars (Krabbe et al. 1995), as well as many less 
massive O and probably B stars (Eckart et al. 1999). The 
Hel stars are believed to be post-main sequence supergiant 
stars of ~ 20 — 100 M Q , close to the Wolf-Rayet stage (Na- 
jarro et al. 1997). The most massive of these stars have a 
total age of < 3Myr, while the less luminous stars could 
have ages up to ~ 8Myr. Krabbe et al. (1995) have ar- 
gued that the most likely origin of these stars is in a small 
starburst ~ 3 — 7Myr ago, in which > 10 4 M Q of young 
stars were created. 

In situ formation of these stars is problematic, however, 
because of the strong tidal field of the Galactic nuclear 
bulge and central black hole. The tidal forces from the 
Galactic nucleus alone are sufficient to unbind gas clouds 
with densities n cr j t < 10 7 cm _3 (1.6 pc/ifo) at galacto- 
centric radii Rq (e.g., Morris 1993), while clouds in the 
nuclear gas disk have densities of 10 4 — 3 x 10 5 cm~ 3 (Gen- 
zel 1989, Glisten 1989). Hence Morris (1993) has argued 
that the formation of the central star cluster must have 
been externally triggered to achieve the required high den- 
sities, perhaps in a cloud collision. Two clouds coming in 
from ~ 10 pc would need to collide near the gravitational 
radius of the black hole (~ 0.5 pc), at a velocity of several 
hundreds of kms -1 ; it is unclear whether the required 
densities and a high efficiency of star formation could be 
reached. Another model in which the massive stars form 
through collisions and mergers of lower mass stars in the 
high-density nuclear cluster now appears unlikely, both be- 
cause too few massive stars would form (Lee 1994), and 
because similar stars have been found in the Arches and 
Quintuplet clusters some 30 pc away from the center (Na- 
gata et al. 1995, Cotera et al. 1996, Figer et al. 1999a). 

This letter therefore explores an alternative idea: that 
the young stars now seen in the Galactic Center (GC) 
formed further out in a massive star cluster that subse- 
quently spiralled into the nucleus and tidally dissolved. 
One of the exciting results from HST has been the dis- 



covery of young star clusters in a variety of starburst en- 
vironments (e.g., Whitmore & Schweizer 1995, O'Conncll 
et al. 1995, Oestlin et al. 1998). We also know now that 
nuclear star clusters are common in spiral galaxies (Car- 
olio et al. 1998, Matthews et al. 1999). The Arches and 
Quintuplet clusters at ~ 30 pc distance from the GC tes- 
tify that a similar star formation mode has been occurring 
in the nuclear disk of the Galaxy. Both clusters have ages 
of a few megayears and estimated total masses (extrapo- 
lating the IMF to 1 Mq) of - 10 4 M (Figer et al. 1999b). 
The orbits of massive clusters evolve by dynamical fric- 
tion against field stars in sufficiently dense stellar systems 
(Tremaine, Ostriker & Spitzer 1975). Here I show that a 
massive cluster formed well outside the central few parsecs 
will indeed spiral into the center within the lifetime of its 
most massive stars, losing much of its mass on the way 
inwards until finally its dense core is disrupted deep in the 
nucleus. 

2. MASSIVE CLUSTER INFALL 

Suppose a massive star cluster is formed at initial galac- 
tocentric radius Ri — 30i?3o pc. The mass distribu- 
tion of the nuclear bulge in 2pc < Rq < 30 pc can 
be approximated as an isothermal sphere (Genzel, Hol- 
lenbach & Townes 1994, Fig. 7.1), with circular velocity 
v c = 130t>i3o kms -1 in this range of radii, and 

(1) 
30 pc 

corresponds to 10 5 ' 5 — 10 4 ' 5 H-atoms per cm 3 , so in this 
region molecular clouds with the densities observed in the 
GC may indeed collapse and form stars. 

The newly formed cluster will lose orbital energy by dy- 
namical friction against the nuclear bulge stars and spiral 
into the center on a time-scale 

r * = MolDv) £ = 3A x ^^WAio 1 yr, (2) 

where the cluster mass is m c — 10 6 mg M Q and Aio = 
ln(0.4iV)/10 > 1 with N the effective number of field 
stars (Binney & Tremaine 1987). Thus a massive clus- 
ter, m e ~ 10 6 M Q , will indeed spiral into the center from 



M G (R G ) = 3.9 x 10 6 ^ 30J R G [pc] M G 
Eq. (1) implies that the mean density within 10 
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Ri < 30 pc within the lifetime of its most massive stars. 
Tdf corresponds to only a few initial rotation times around 
the center; at Ri 

r rot = 2ttRJv c = 1.4 x 10 6 i? 3 owr 3 o y r - (3) 
The friction time-scale Tdf is dominated by the time 
spent at large radii. In these initial phases the cluster 
will be embedded within its parent molecular cloud while 
both are dragged inwards together. Because the part of the 
molecular cloud not converted into stars is conceivably the 
major part of the total initial mass, the dynamical friction 
on the cloud (Stark et al. 1991) contributes to the orbital 
evolution of the cluster. Tdf is shorter or comparable to 
the time-scales of other processes that might separate the 
surrounding gas cloud from the cluster. Magnetic friction 
(Elmegrecn 1981) brakes the gas cloud in ~ 3 x 10 7 yr 
(Morris 1996), assuming a mG field continuously perme- 
ates the cloud and rotates more slowly than the cloud by 
a fraction of the circular velocity. Mass loss from bulge 
stars (Jenkins & Binney 1994), even if it can cool and set- 
tle into the disk plane, amounts to only ~5x 10 4 M over 
3 x 10 6 yr into the central 30 pc, so dilutes the cloud angu- 
lar momentum only slightly. Most dangerous for the cloud 
is the radiation from the massive stars, which drives an ex- 
panding HII region; e.g., for the 6 x 10 51 Ly c photons per 
second from a 3 x 10 5 M cluster (Efstathiou et al. 2000), 
and a cloud of 1.7 x 1O 6 M with density 10 4 - 5 cm _3 , the 
HII region front will reach the surface after 2.5 x 10 6 yr 
in standard spherical theory (Spitzcr 1978). Thereafter 
the gas expands and is tidally stripped, but will take a 
further > 1-2 orbital times before it also gravitationally 
decouples from the cluster. To separate most of the cloud 
from the cluster will take longer if the cluster is located 
near the cloud surface or if the cloud is fractal-like. Thus 
the surrounding molecular cloud could contribute to the 
dynamical friction for ~ Tdf, certainly in the crucial initial 
phases, and is likely to be eventually left behind. In detail 
this depends on the parameters of the cloud and cluster, 
but the most massive clouds are the most robust. 

As it spirals to the Galactic center, the parent cloud and 
embedded cluster will constantly lose mass because of the 
strong tidal field. We can estimate the effect of this mass 
loss on the friction time with a simple model in which the 
initial mass m ci of the combined cloud and cluster is dis- 
tributed according to an isothermal profile, tidally limited 
at radius r = ru- 

m c (r) = m ci (r/r ti ), (4) 

Tti = ( jSk^ ) Rl = 6 - 2m 6 /3w i30 /3i? 30 3 P c - (5) 

For comparison, the half-mass radius of the present Arches 
cluster (~ 1O 4 M ) is - 0.2 pc (Figer et al. 1999a). As the 
cluster spirals in, successive outer shells of the mass dis- 
tribution are peeled off by the tidal field. The division be- 
tween cluster and cloud need not be specified now, but it is 
clear that at some stage the cloud envelope will have been 
removed and the cluster proper will begin to be stripped. 
If Rq is the current galactocentric radius, the tidal radius 
decreases according to 

r t /R G = ru/Ri, (6) 
and if internal evolution is neglected, the cluster mass de- 
creases as 

m c {R G ) = {R G /Ri)m ci . (7) 



When the entire cluster is finally disrupted at radius Rdis, 
only a fraction <~ R^/Ri of the initial mass m C i will have 
arrived at i?dis- 

We can now write a modified dynamical friction equa- 
tion using the same arguments as in Binney & Trcmaine 
(1987), and assuming that the instantaneous specific 
angular momentum loss due to the frictional force is 
the same for stripped material and material that re- 
mains bound to the cluster. This gives RcdRc/dt = 
-0A28Gln(0AN)v^ 1 m cl (R G /R l ), the solution of which 
is 

r df (Ri, R G ) = R l R G v c /0A28 In(0.4JV) Gm a (8) 
= 6.2 x 10 6 (i? G /i? J )^30«i3om 6 - 1 A 10 1 yr. 

Thus in this simple model the total time to spiral in 
from radius Ri is just twice that when the mass m ci re- 
mains constant [eq. (2)], and the time taken from radius 
Rq < Ri for a cluster that started out tidally limited at 
Ri is linearly proportional to Rq. To reach the center in 
~3x 10 6 yr from Ri ~ 30 pc (10 pc), a mass-losing cloud- 
cluster system must have an initial mass m ci ~ 2 x 10 6 M 
(~ 2 x 10 M©). The core of a cluster with initial mass 
2 x 10 6 M formed at Ri ~ 10 pc would have spiralled into 
the center after 300,000 years. These time-scales may be 
slightly overestimated because we have neglected torques 
from previously stripped material and its wake, left be- 
hind on orbits of lower frequency, as well as any additional 
drag from the nuclear gas disk. Perturbed material within 
a factor 1.5-2 in galactocentric radius contributes to the 
frictional drag (Tremaine & Weinberg 1984); over this ra- 
dial range the cluster loses about half of its mass [eq. (7)]. 

One constraint is that the cluster must not evaporate in 
the strong tidal field of the nuclear bulge before it reaches 
the center. In the m c = 2 x 1O 4 M evolutionary models 
for the Arches and Quintuplet clusters by Kim et al. (1999) 
the evaporation time is several Myr. Evaporation occurs 
on a time-scale proportional to the half-mass relaxation 
time 

t rh = 2.0 x 10 8 m\ /2 rf 2 m" A^ 1 yr, (9) 

(Spitzer & Hart 1971) but depends also strongly on the 
strength of the tidal field. Here = n pc is the clus- 
ter half-mass radius and = m* jQ M© the average stel- 
lar mass. t r h is sensitive to the uncertain half-mass ra- 
dius. For the Arches cluster Figer et al. (1999a) esti- 
mated rh — 0.2 pc from the observed high-mass stars, i.e., 
without correction for initial or dynamical mass segrega- 
tion. We can rewrite t r h oc rnj 2 ^ 2 cx to c( o~ 1//2 and con- 
sider the mean density to be fixed by the external tidal 
field. This suggests independently that only clusters sig- 
nificantly more massive or formed at smaller radii than the 
Arches cluster can reach the central parsec. 

3. FINAL TIDAL DISRUPTION AND THE HE I STAR 
CLUSTER 

The preceding discussion shows that only the inner parts 
of a massive young star cluster will spiral into the Galactic 
center proper, while its outer parts and the cloud enve- 
lope will be left behind and distributed further out. There 
are good reasons to believe that the massive stars will be 
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concentrated towards the cluster core. There is some evi- 
dence that high-mass stars form predominantly in or near 
the cores of young clusters (Fischer et al. 1998, Bonnell & 
Davies 1998). After star formation is complete, dynam- 
ical mass segregation will further accentuate the central 
concentration of massive stars, acting on their two-body 
relaxation time-scale. Thus the most massive cluster stars 
end up in the GC, as is needed to explain the observed 
distribution of Hel stars. 

Under what conditions is the density of the cluster core 
high enough to reach the central parsec? The GC Hel 
stars are observed at galactocentric radii Rq = 1 — 10" = 
0.04 — 0.4 pc and their proper motions and radial velocities 
show a coherent rotation pattern (Genzel et al. 2000). In 
the present scenario this would be interpreted as tracing 
the last orbit of the cluster core before final disruption, 
giving a radius of disruption i? d ; s 0.4 pc. In this radial 
range the tidal force of the central black hole dominates 
that from the nuclear bulge. We can thus estimate the 
required mean density of the cluster core at i? d j s as 

p dis = 3m c (i? dis )/47rr t 3 (i? dis ) = 6M./4nR 3 dis (10) 

= 2.2 x 10 7 M 3 (i? dis /0.4pc) -3 M pc~ 3 , 

where we have written the black hole mass as M, = 
3 x 10 6 M 3 M Q (Genzel et al. 2000). Notice the sensitivity 
to the disruption radius. 

The observed average stellar density in the Arches 
cluster, extrapolated to include stars down to 1M Q , is 
6.3 x 10 5 M Q /pc 3 (Figer et al. 1999a). The 2 x 1O 4 M 
models of Kim et al. (1999) show that dynamical evolu- 
tion and mass segregation lead to the formation of a core 
of high mass stars with density reaching p c ~ 10 7 M Q / pc 3 
in the mild core collapse stage, after ~ 1 Myr evolution. 
Central densities of ~ 10 7 M Q pc~ 3 can thus be reached 
in later evolutionary stages, high enough for the core to 
survive to i? d is ~ 0.4 pc, but significantly larger densities 
(required for smaller i? d i s ) would seem problematic. 

The time-scale for this evolution is a small number of 
relaxation times for the massive stars (Chernoff & Wein- 
berg 1990). If evolution is too fast, the core may reexpand 
before the cluster has time to reach the GC; this would 
be the case for the Arches cluster models. For clusters 
with initial masses m c > 3 x 10 5 M Q , the mass segregation 
time for 30 M Q stars is > 2 x 10 6 yr, using (9) and scaling 
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radii such that cx m c ■ Thus in the initial stages dy- 
namical evolution is relatively slow. As the cluster spirals 
in and loses mass, however, the evolution accelerates: the 
half-mass relaxation time decreases because both the clus- 
ter mass and half-mass radius decrease, and the external 
tidal field becomes stronger, so that the mass evaporat- 
ing per relaxation time increases. Thus the most rapid 
dynamical evolution is expected to take place when the 
stripped-down cluster approaches the center. At the time 
of final disruption, its mass is a small fraction of the ini- 
tial mass, comparable to the model clusters studied by 
Kim et al. (1999) in which the time-scale for core collapse 
is ~ 10 6 yr. Thus it is not improbable that such a clus- 
ter reaches core collapse shortly before it arrives at the 
Galactic Center. The high cluster core density required 
for the core to reach i? d j s ~ 0.4 pc would then most eas- 
ily be maintained in a phase after core collapse when the 
energy loss of the core from collisions is compensated by 



the energy gain from three-body binaries. This requires 
the core to be dominated by the 100 most massive stars 
(Binney & Tremaine 1987), previously concentrated into 
this volume by mass segregation. 

After the high-density cluster core has reached i? d ; s ~ 
0.4 pc and begins to disintegrate, it continues to spiral in- 
wards until torques from its friction wake and the material 
lost previously (see §2) become ineffective, that is, until 
the debris has spread in angle by A<p ~ ir. The time for 
the debris within ar t of the cluster center to spread by 
A(f) = tt is 



^spr — 



R(li; 



dfi 



di? G 



^di 



2ar t (R dis ) 



trot 

6 ar t (Rdis) 



(11) 



which gives of order t rot /a according to eqs. (5), (6). 
The rotation period near the black hole is t rot = 
13700M 3 " 1/2 (i? G /0.4pc) 3 / 2 yr, so for the cluster center 
t Iot /a can be several 10 4 yr. For comparison, the dynami- 
cal friction time (2) for 10 4 M Q to spiral from i? d i s = 0.4 pc 
into the center completely is only 50000 yr. This suggests 
that the debris of the cluster core will spiral inwards signif- 
icantly even after the final disruption, and that therefore 
the rotation pattern observed down to Rq — 0.15 pc is 
consistent with disruption at radii ~ 0.4 pc. It will be 
interesting to simulate this disruption event in the black 
hole's tidal field with numerical star cluster models. 

4. DISCUSSION 

In summary, the previous sections have shown that un- 
der plausible conditions a young star cluster formed in the 
Galactic nuclear disk may spiral into the GC within the 
life-time of its most massive stars, and that the density in 
its core of high-mass stars can be such that it would ac- 
tually reach and dissolve in the region where the GC Hel 
stars are observed. The important parameters are the ini- 
tial galactocentric radius, the cluster mass and radius, and 
the mass and density of the surrounding molecular cloud. 
To reach the center from initial radii Ri ~ 10 — 30 pc, 
the cluster must be substantially more massive than the 
Arches and Quintuplet clusters or the frictional force must 
at least initially be augmented by the wake of the sur- 
rounding molecular cloud, or both. Large cluster masses 
are also required to avoid too rapid dynamical evaporation 
and evolution in the strong tidal field. 

A prediction of the present model for the origin of the 
GC Hel stars is that the massive stars in the core of the 
disrupting young cluster should end up at smaller galacto- 
centric radii than lower mass stars in the cluster envelope. 
Krabbe et al. (1995) estimate the total stellar mass corre- 
sponding to the nuclear Hel star cluster as ~ 1.5 x 10 4 M , 
assuming a mass function N(m) dm oc m~ 2 dm down to 
m* = 1 M Q . The number of low-mass stars associated 
with the cluster is unknown, but a standard Salpeter IMF 
with the same number of the highest mass stars, extrap- 
olated to m* = 0.1 M Q , would give ~ 10 5 M Q . Most of 
these stars should now be found at galactocentric radii of 
a few parsec (some tens of arcsec) or further out if the 
cluster was even more massive. Could these stars be re- 
lated to the young population inferred by Philipp et al. 
(1999), or was the formation of the cluster just part of a 
larger-scale starburst? 
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As in alternative models, the observed counterrotation 
of the Hel stars with respect to Galactic rotation (Gen- 
zel et al. 2000) is not inherent in the present scenario. 
However, a large molecular cloud entering the nuclear disk 
with orbital angular momentum roughly antialigned with 
Galactic rotation would sooner or later suffer a strong colli- 
sion, and this in fact might be the event triggering the for- 
mation of the cluster. The subsequent dynamical friction 
on the cluster and its surrounding cloud envelope would 
then tend to circularize the counterrotating orbit while the 
cluster spirals in, explaining the rotating torus structure 
of the Hel stars observed now. This observation would be 
more difficult to explain in models where the GC young 
stars form in a cloud collision near their present location, 
because some of these stars should then still reflect the 
cloud's initial orbital motion. 

How common could cluster accretion events be in the 
GC? The presence of ~ 10 8 yr old AGB stars in the central 
parsec (Krabbe et al. 1995) suggests that it has happened 
before, and perhaps one of the massive molecular clouds 
in the vicinity is the next candidate. Clearly, a cluster 
accretion rate of ~ 10 5 M Q per 3 x 10 6 yr would build up 
substantial mass in the nuclear bulge over time. The total 
mass inflow rate of ~ 0.1 Mq/ yr inferred from dynamical 
friction on the GC giant molecular clouds on a^-orbits at 
~ 100 pc (Stark et al. 1991), and somewhat more from gas 
flow through the ILR at ~ 200 pc (Gerhard 1992, Morris 
& Serabyn 1996) would supply sufficient gaseous material 



to maintain this rate of cluster formation in addition to 
other distributed star formation. Detailed observations 
of the young stellar population in the GC are needed to 
see whether the process actually happens recurrently and 
whether it is a substantial factor in building up the nuclear 
bulge. Early-on, clusters would spiral in as far as there are 
field stars to provide the background for dynamical fric- 
tion. Later clusters would come in further and build up 
a density gradient (assuming the angular momentum can 
be further transported outwards), but only some would 
reach the vicinity of the black hole. In the simple model 
described in Section 2, where both the nuclear bulge and 
cluster density profiles are oc r~ 2 , the ingoing cluster's 
tidally stripped material would also be distributed oc r~ 2 . 

Independent of the rate of cluster accretion, the recent 
arrival in the GC of a cluster with a large number of mas- 
sive stars would have profoundly changed the physical con- 
ditions in the central parsec and may in fact be responsible 
for the present lack of activity of Sgr A*. If so, stellar pop- 
ulation studies of the frequency of cluster accretion in the 
GC will have some general impact on understanding the 
nuclear activity cycles in spiral galaxies. 
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ported by grant 20-56888.99 from the Swiss National Sci- 
ence Foundation. 
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